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ABSTRACT 

We present K-band imaging of two ~ 30' x 30' fields covered by the MUSYC 
Wide NIR Survey. The 1030 and 1255 fields were imaged with ISPI on the 
4m Blanco telescope at CTIO to a 5a point-source limiting depth of K ~ 20 
(Vega). Combining this data with the MUSYC Optical UBVRIz imaging, we 
created multi-band K-selected source catalogs for both fields. These catalogs, 
together with the MUSYC K-band catalog of the ECDF-S field, were used to 
select K < 20 BzK galaxies over an area of 0.71 deg 2 . This is the largest area 
ever surveyed for BzK galaxies. We present number counts, redshift distributions 
and stellar masses for our sample of 3261 BzK galaxies (2502 star-forming (sBzK) 
and 759 passively evolving (pBzK)), as well as reddening and star formation rate 
estimates for the star-forming BzK systems. We also present 2-point angular 
correlation functions and spatial correlation lengths for both sBzK and pBzK 
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galaxies and show that previous estimates of the correlation function of these 
galaxies were affected by cosmic variance due to the small areas surveyed. We 
have measured correlation lengths r of 8.89 ±2.03 Mpc and 10.82 ± 1.72 Mpc for 
sBzK and pBzK galaxies respectively. This is the first reported measurement of 
the spatial correlation function of passive BzK galaxies. In the ACDM scenario 
of galaxy formation, these correlation lengths at z ~ 2 translate into minimum 
masses of ~ 4 x 1O 12 M and ~ 9 x 1O 12 M for the dark matter (DM) halos 
hosting sBzK and pBzK galaxies respectively. The clustering properties of the 
galaxies in our sample are consistent with them being the descendants of bright 
LBG at z ~ 3, and the progenitors of present-day > 1L* galaxies. 

Subject headings: galaxies: high-redshift, evolution, statistics - - cosmology: 
large-scale structure of universe — catalogs — surveys 



Introduction 



The processes involved in the formation of galaxies and the phenomena that drive their 
consequent evolution are among the most important and fundamental problems in astron- 
omy. In order to understand these processes we need to study the properties of statistically 
significant samples of galaxies throughout cosmic time. During the last decade there have 
been enormous advances in this direction. The advent of 8m class telescopes, a new genera- 
tion of wide field optical imagers and infrared detectors, and a series of different photometric 
selection techniques, have provided us with the ability to select large samples of interme- 
diate and high redshift galaxies. Probably the most successful method applied to date in 
order t o build these samp les has been the Lyman Break Galaxy (LBG) dropout technique at 
z ~ 3 (jSteidel et al.lll996l ) which has recently been extende d to lower redshifts (z ~ 2) using 



U„.G1 Z two-color photometry, the BM/BX galaxy selection (Erb et al.ll2003l ; lAdelberger et al. 



20041 ) . Othe r types of z > 1 galaxies include narrow - band selected Lyman A lpha Emitters 



(LAE) (e.g. iHu et all (119981 ): iGawiser et all fl2006ri ): iGronwall et all (12007^. Distant Red 



Galaxies (DRG) selected by having red near-infrared colors J — K > 2.3 (IFranx et al.l 12003 



van Dokkum et al.l 120031 ). Extremely R ed Objects (ERO) selected by their red op tical to 



near-infrared colors, usually R — K > 5 ([Thompson et al.lll999t iCimatti et al.ll2002al) as we ll 



as far-IR and sub- mm detected, highly obscured galaxies at high redshift (ISmail et al.ll2000l ). 



Another selec tion technique, of great importance to this work, was introduced by 
Daddi et al.l (120041 ) who defined a two color criteria based on BzK photometry of K-selected 
galaxies, which is capable of identifying both passive (pBzK) and star-forming (sBzK) galax- 
ies in the 1.4 < z < 2.5 range, and to distinguish between the two populations. Among the 
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available NIR bands, the K band is the best choice to detect intermediate and high redshift 
galaxies from the ground. It is the reddest band for which sensitivities and angular reso- 
lutions comparable to the ones obtained by ground-based optical imaging can be achieved. 
Since galaxies can be selected efficiently by their rest-frame optical emission redshifted to the 
NIR bands, K-selection gives us access to the "spectroscopic redshift desert" (1.4 < z < 2) 
where galaxies show no distinctive spectral features in their observed optical spectra. One 
of the main advantages of the BzK technique is that it is almost insensitive to the amount 
of reddening in the star-forming galaxies. NIR observations allow us to detect and iden- 
tify distant galaxies in a more unbiased way than optical imaging, since the effects of dust 
extinction are lower at these wavelengths and the observed light is dominated by evolved 
stellar populations, making detection more independent of spectral type. NIR luminosity is 
directly related to the stellar mass of the galaxy. The BzK selection technique provides us 
with a powerful tool to explore the z ~ 2 universe in a relatively unbiased way. 

Despite all these advances, and the fact that we have converged into a single cosmo- 
logical model, the assembling and evolution of galaxies remains poorly understood. We still 
lack a comprehensive picture of galaxy populations at intermediate and high redshifts. All 
the selection techniques mentioned above only give partial views of the universe at specific 
epochs, and the biases involved in their selectio n, as well as the overlapping between differ- 



ent techniques is jus t starting to be understood (IReddv et al.ll2005l ; Ivan Dokkum et al.ll2006 



Quadri et al.ll2007bl ; lLane et al.l 120071 ; iGrazian et al.l 120071 ). In particular, optically selected 
samples require galaxies to be bright in their rest-frame UV, and are therefore biased against 
obscured systems in which the UV luminosity of a star-forming galaxy is quenched by the 
dust content of its ISM, as well as against passively evolving systems with little intrinsic 
UV emission. Hence, studies based on these samples constrain the behavior and evolution 
of unobscured star formation activity, more than the actual stellar mass assembly of these 
systems. At z ~ 2 the most popular select ion techniques used so far include UV selected 
BM/BX galaxies, DRGs, and BzK galaxies. IReddv et all J2005h showed that BM/BX selec- 
tion misses 40% of star-forming BzK galaxies whe n both sample s are su bject to a common 
NIR magnitude limit of K < 20. Recent work by I Quadri et al.l (l2007bl ) estimated a larger 
missing fraction o f 60% to K < 20, w hich decreases to a 35% when going to deeper NIR 
fluxes (K < 21). IGrazian et al.l (120071 ) computed an even larger fraction of missed sBzK 



galaxies of a 62% at a deeper NIR magnitude limit of K < 22. In all the mentioned works 
BM/BX galaxies completely miss passive BzK galaxies at z ~ 2. UV based s election tech- 
niques perform even more poorly when trying to recover the DRG population. iQuadri et al. 



fl2007ri ) showed that only 32% of DRG in the K < 21 MUSYC Deep NIR sample have colors 
consistent with the BM/BX/LBG criteria, and only 14% also fulfill the TZab < 25.5 cut used 
in most ground-based optically selected samples. This is in agreement with the 12% recovery 
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fraction estimated bv lReddv et al.l (120051) to the sa me NIR magnitude limit and higher than 
the 3% fraction computed by iGrazian et al. (bo07h for a deeper K < 22 DRG sample. 



On the other hand, the sBzK criteria successfully recovers 80% o f the z > 1.4, K < 20 
BM/BX galaxies in Reddy's sample, in agreement with the results of iQuadri et al.l (j2007bl ) 
for a K < 2 1 sample. This fractio n improves to up to a 99% when going to a fainter NIR limit 
of K < 22 (iGrazian et al.l 120071 ). The BzK technique is also much better than rest-frame 
UV se lection techniques at recovering the DRG population. For a bright K < 19.3 sample of 
DRG, lLane et al.l (120071 ) obtained a recovery fractio n of 55% being selec ted as either sBzK 
or pBzK. A higher fraction of 65% was obtained by iReddv et al.l (120051 ) for K < 21 DRG, 



again in agree ment whith the result s of IQuadri et al 



(120070 ). When going to fainter NIR 



fluxes(-K" < 22) IGrazian et al.l (120071 ) estimated an almost complete recovery fraction of 99%. 
These numbers show that the BzK technique is the most inclusive galaxy selection method 
at z ~ 2, mostly because is able to include galaxies that span a wide range of both optical 
and optical-to-NIR colors. 



Previous samples of BzK galaxies include the seminal work by iDaddi et al.l (12004 ) who 



selected K < 20 BzK galaxies over an area of 52 arcmin 2 , a survey by iKong et al.l (120061 ) 
who covered 320 arc min 2 to K < 20 a nd 920 arcmin 2 to a shallower depth of K = 19, and 
more rece nt work by lLane et al.l (120071 ) who surveyed 0.56 deg 2 to K = 20.6 in the UKIDSS 
field, and lHayashi et al.l (120071 ) who selected very faint BzK galaxies (K = 21.3) over a 180 
arcmin 2 area in the Subaru Deep Field. 

In this work, we present the MUSYC Wide NIR Survey, and the clustering and physical 
properties of BzK galaxies selected over an area of 0.71 deg 2 to a limiting depth of K = 20. 
This is the largest area ever surveyed for BzK galaxies, and the largest sample of these objects 
available to this depth. In Section 2 we introduce the MUSYC survey and present the Wide 
NIR Survey design. Sections 3 and 4 discuss the observations and data reduction. Source 
detection, photometry and the procedures used to build the K-selected source catalogs are 
discussed in Section 5. In Section 6 we explain our Star-Galaxy classification method, and 
we present K-band selected galaxy number counts in Section 7. We present the details of 
the BzK selection technique used in this paper in Section 8. The results coming from the 
analysis of the selected BzK sample, redshift distributions, physical properties and clustering 
are presented in Section 9. Finally, we summarize our results and present conclusions in 
Section 10. 

All optical magnitudes are given in the AB system, and NIR magnitudes in the Vega 
system. Throughout the whole paper we adopted the standard cosmological parameters, 
H Q = 71 km s _1 Mpc _1 and Qm = 0.27, = 0.73 and as = 0.8. All quantities reported are 
comoving, so correlation lengths scale as h? , number densities as hj and halo masses as 
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hjQ, unless explicitly stated. 



2. The MUSYC Wide NIR Survey 



The MUltiwavelength Survey by Yale Chile (MUSYC) consists of deep optical (UBVRIz) 
and near-infr ared imaging o f four 3 0' x 30' fields. The survey design and optical images are 
described by iGawiser et al.l (j2006af ). The four fields (Extended Hubble Deep Field South 
(EHDF-S), Extended Chandra Deep Field South (ECDF-S), SDSS1030+05 and Castl255) 
are at high galactic latitudes (| b |> 30) and show low 100 /iin emission and Nh in order to 
minimize contamination by stars and the Galaxy ISM. 

The near-infrared imaging of the MUSYC fields forms part of two complementary cam- 
paigns. The MUSYC Deep NIR Survey imaged four 10' x 10' sub-fields in three NIR bands 
to typical 5a limiting d epths of J ~ 2 2.5, H ~ 21.5 and K ~ 21. The survey design and 
results are presented by iQuadri et al.l (j2007bl ). The MUSYC Wide NIR Survey consists of 
K-band imaging of three full 30' x 30' MUSYC fields (ECDS-S, SDSS1030+05 and Castl255) 
to a limiting depth of K ~ 20 as well as J and H-band imaging of the ECDF-S. 

In this work we present the NIR imaging and K-band selected source catalogs for two 
fields of the Wide NIR Survey: SDSS 1030+05 and Cast 1255 (hereafter 1030 and 1255). 
Table [1] gives relevant informa tion about th e surveye d fields. The wide NIR imaging of the 



ECDF-S will be presented by iTaylor et al.l (Iin prep.1 ). H-band imagin g to similar depth of 



EHDF-S is already available from the Las Campanas Infrared Survey (jChen et al.l 120021 ) . 
MUSYC also includes narrow- band 5000A imaging of the whole s quare degree to search 



for Lyman Alpha E mitters and Lyman Alpha Blobs at z — 3.1 (IGawiser et al.l 12006b 



Blanc et a 



m prep. 



), mid-infrar ed from the Spitzer IRAC MUSYC Public L egacy Survey 



(SIMPLE. lLabbe et al.l (lin prep.l )). Chandra X-ray imaging (jVirani et al.ll2006l ) of the ECDF- 
S, XMM-Newton X -ray observations of 1030, and spectroscopic follow up of the four fields 
(ILira et al.l lin prep.l ). 



More information about MUSYC, together with publicly available data can be found at 



http:/ '/www. astro, yale.edu/musyc/ 



3. Observations 



Our K-band images of the 1030 and 1255 fields were taken on the nights of 2003 January 
31, 2003 May 7-12, 2004 February 11-14, 2004 April 10-13, 2005 April 14-22, and 2006 April 
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18-23, using the Infrared Sideport Imager (ISPI) on the 4m Blanco Telescope at CTIO. The 
detector is a 2Kx2K HgCdTe HAWAII-2 array which delivers images with a 0.305" pixel -1 
scale over a field of view (FOV) of 10.5' x 10.5'. 

Since the FOV of ISPI is smaller than the 30' x 30' area covered by MUSYC optical 
images, we divided each field in nine sub-fields (NE, E, SE, S, SW, W, NW, N, C) that we 
later mosaic as described in Section 4.4. 

The dominance of the sky background and its rapid variability determines the observing 
strategy in the K band. Dithered short exposures are necessary in order to subtract the 
background and keep the number of counts in the linear regime of the detector. Individual 
raw frames are the result of co-adding short exposures. Our science frames have integration 
times of 4 x 15, 5 x 12 or 6 x 10 seconds, depending on the intensity of the sky background 
at the moment of the observations. We used a quasi-random dither pattern inside a 45" box. 

Dome flats and dark frames were obtained at the beginning and the end of every night. 
Flats are built by taking the difference between the combined flat frames taken with the 
lamp on and the combined frames taken with the lamp off in order to subtract the thermal 
background. 

The K' filter on ISPI was replaced in April 2004 by a new Ks filter. Around 20% of our 
observations were completed using the K' filter, and the rest of the survey was completed 
using the new Ks filter. Figure Q] shows the transmission curves of the two filters after 
multiplying by the detector quantum efficiency and the atmospheric transmission at CTIO. 
The differenc e in effective waveleng th for the two filters is ~ 30A, negligible compared to the 



filters width. lOuadri et al.l (l2007bl ) calculated a conversion between both filters of 0.02 mag, 
smaller than our typical photometric uncertainties. Therefore in the following discussion we 
do not distinguish between the two filters and we simply refer to them as the K filter. 



4. Data Reduction and Image Properties 

4.1. Image Processing 

Data reduction was performed using standard IRAe3 routines, modified IRAF rou- 
tines and custom IDL tasks. Our reduction scheme is based on the methods described by 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomylnc, under cooperative agreement of the National Science 
Foundation. 
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Labbe et all (120031 ) and lQuadri et al.l fl2007bf ). The external IRAF package XDIMSUIvflwas 



used throughout most of the steps of the reduction. 

We subtracted dark frames of the same exposure time and number of co-adds from our 
science images. Flat frames were constructed for each night as described above, but flat 
fielding was not performed until after background subtraction. This approach is adopted 
because the multiplicative nature of the flat fielding process scales up the contribution to 
the noi se coming fro m the uncertainty in the flat frame by the number of counts in the science 
image ( Uoyce!ll992l ). This can be shown by deriving the uncertainty in a given pixel after 
background subtraction and flat fielding. We consider the two cases where (1) background 
subtraction is performed before flat fielding and (2) when flat fielding is performed before 
background subtraction. 

Consider a pixel p which has a background flux b and object flux o in the image un- 
der construction, but only background flux in the other N frames used to estimate the 
background. In case 1, the value of the pixel after background subtraction is 



and the uncertainty is 



p = b + o— <b> 



Jl jz , ji , a b 
a p = a b +a Q + — 



(2) 



now, we divide by the flat and we get 



V 



b + o- <b> 
1 



(3) 



and the final uncertainty in case 1 becomes 



a, 



a\ + al + J- (b + of a) _ 2{b + o) <b> o) < 6 > 2 a) 

PA n A 



f 2 f 4 f 4 

Now, in case 2, we first divide by the flat obtaining 



.Z 1 



(4) 



2 XDIMSUM is the Experimental Deep Infrared Mosaicing Software package developed by L. Davis, and 
it is available at http://iraf.noao.edu/iraf/ftp/extern/xdimsum 
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b + o 



(5) 



P = 



f 



with an uncertainty 




+ 



{b + ofa} 
T 4 



(6) 



and then subtract the background obtaining 



P = 



(b + o) <b> 

~1 T 



(7) 



with a final uncertainty for case 2 of 




+ 



ib + ofa) al 



+ 




(8) 



/4 Nf 2 



It can be seen that the noise in case 2 is higher than in case 1 by an amount 2(6 + 6) < 
b> a'j/f 4 . Since the typical background level in our images is b ~ 30000 counts, and after 
sky subtraction is ~ counts, it is clear that the contribution from the uncertainty in the 
flat frame will be strongly suppressed by flat fielding after sky subtraction. In the above 
derivation we have neglected the contributions of dark current and read noise. 

Background subtraction was performed using a two-phase process that avoids overes- 
timating the background around bright objects. A background image is created for every 
science frame by combining eight dithered frames adjacent in time. The combination is an 
average excluding the lowest and highest value at every pixel position. These background 
frames are subtracted from the science images to create preliminary background subtracted 
images. The positions of six stars are used to calculate the relative offsets between dithered 
background subtracted frames and the images are then combined after registering them to 
a common reference using fractional pixel shifts. 

Objects are then detected in the combined image using a thresholding algorithm and an 
object mask is created. The mask is deregistered to the reference of each original dithered 
science frame and the background subtraction process is repeated. In this second phase, 
the background image for each science frame is created as above but masking the detected 
objects in the time adjacent frames before combining them. This background frame is then 
subtracted from the science image. 
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A second set of stars in these final sky subtracted images is used to measure the FWHM 
and calculate the relative flux scale between the images. The central 300 x 300 pixel region 
of each frame is used to calculate the pixel-to-pixel rms of the background. These values are 
used to compute weights used to optimize the signal-to-noise in the seeing disk during the 
combination of the background subtracted images. The weights are given by, 



* (flux_scalei ■ rmsi ■ FWHMi) 2 

A bad pixel mask (BPM) is created by flagging all 4cr deviant pixels in the flat. XDIM- 
SUM uses a cosmic ray detection algorithm and combines this information with the BPM 
to create a rejection mask for each science frame. Satellite trails found during the visual 
inspection of the images are added by hand to the corresponding rejection mask. The ISPI 
detector is not able to properly flush extremely bright pixels during readout, retaining mem- 
ory of previous exposures in the cores of bright objects. To avoid residuals from previous 
exposures we construct a previous frame bright object mask for each science image. 

All background subtracted frames were visually inspected in order to reject bad quality 
images. There are two main reasons for rejection. The first one is catastrophic background 
subtraction on images due to extreme variations of the background. The second cause of 
rejection is a poor shape of the point spread function (PSF) due to focusing and tracking 
problems. Around 10% of our images were rejected because of these problems. All the 
accepted background subtracted images are then flat fielded dividing by the normalized 
dome flat. 

Flat fielded images usually show residuals caused by variability of the bias structure 
of ISPI's detector. We remove these features by fitting a low order polynomial surface to 
each quadrant of the images after masking the objects and subtracting this surface from the 
images. Finally all the images are registered and combined using the above weights, rejecting 
bad pixels, cosmic rays and residuals from previous frames. An exposure time map is also 
generated. Details about the K-band images obtained for each sub-field in 1030 and 1255 
are given in Table [2J 



4.2. Astrometric Projection, PSF Matching and Photometric Calibration 

In order to measure accurate colors, we need to perform aperture photometry in different 
bands over the same physical region projected on the sky. The MUSYC optical images have 
been re-sampled to a tangent plane projection with a uniform pixel scale and show an rms 
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astrometric error of 0.2" ( IGawiser et al.ll2006al ). We re-projected the K-band images so they 



share the same projection, image size and world coordinate system of the optical images. 

We used the standard IRAF tasks GEOMAP and GEOTRAN to compute and apply 
a transformation between the logical coordinate system of the ISPI images and the logical 
system of the reference MUSYC BVR image. For each sub-field image we used the positions 
of ~ 100 stars detected in both the K-band and the BVR image to fit the coefficients of the 
2D 6th order Chebyshev polynomials used for the transformation. The coordinates of the 
stars in the transformed K-band images show an average rms of 0.05" with respect to the 
coordinates in the BVR image, so our astrometric error is dominated by the uncertainty in 
the optical images. After re-sampling the pixel scale changes from 0.305" pixel -1 to 0.267" 
pixel -1 . 

The measurement of accurate colors requires matching the PSF of the optical and K- 
band images. We built an empirical PSF for each K-band sub-field image by registering and 
combining the images of ~40 stars. We used the same method with ~ 120 stars to build 
empirical PSF for the optical images in each band for the two fields. Tables H] and [3] present 
the FWHM of the empirical PSF of the K-band and optical images respectively. Since the 
PSF shows broader wings than a Gaussian profile of the same FWHM, convolution with 
a purely Gaussian kernel is not able to match both the core and the wings of the PSF at 
the same time. This implies the need to use a non-Gaussian kernel to obtain accurate PSF 
matching. We matched the PSF of all the images to a Moffat profile with a FWHM of 1.2" 
for 1030 and 1.3" for 1255. These target PSF were chosen to be slightly broader (both in 
the core and the wings) than the PSF of the image with the largest FWHM in the field. 

We used the IRAF task LUCY to compute the appropri ate convolution kernels for th e 



PSF matching using the Lucy- Richardson iterative algorithm (jRichardsonlll972i ; lLucylll974l ). 
LUCY deconvolves the target PSF into the input PSF and a convolution kernel. The kernels 
are then normalized and convolved with the images to match their PSF to the target Moffat 
profile. In order to quantify the quality of the method we constructed growth curves for all 
the stars used to build the empirical PSF and median combined them to create a master 
growth curve for each PSF matched image. After convolution the master growth curve 
matches the Moffat profile growth curve to better than 1% at all radii for all the images. 



We flux calibrated our images using the 2MASS survey (jSkrutskie et al.ll2006l ). As can 
be seen in Figured] the 2MASS K-band is very similar to the ISPI K-band used in this work. 
We have computed and applied a correction of 0.03 magnitudes to the zero-points to account 
for the difference in the filters, which is of the order of our photometric errors. In each of 
our K-band sub-fields we performed photometry in a 14" diameter aperture of an average 
of 20 stars with 2MASS counterparts and not affected by saturation, non-linearity or close 
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neighbors in the MUSYC imaging. We calculated the zero-points of each sub-field image by 
comparing our 14" aperture magnitudes to the aperture-corrected K-band magnitudes of the 
same stars in the 2MASS point-source catalog (which were measured in 4" apertures and 
corrected using a curve-of-growth measured out to 14"). Uncertainties in the zero-points 
were computed based on the deviation of the zero-point values obtained for single stars and 
are typically ~ 0.02 mag. Table [2] presents the zero-point values for each sub-field in 1030 
and 1255. All sub-fields were scaled to the zero-point with the lowest uncertainty in the field 
so the zero-points of our final K-band images are 22.03 for 1030 and 22.38 for the 1255 field. 



4.3. Correlated Noise Properties 



Derivations of photometric uncertainties typically assume that the background noise 
properties are characterized by Poisson statistics so that the background noise in a given 
aperture is simply given by ^n P i X a P i X , with a P i X the pixel-to-pixel rms and n P i X the number 
of pixels in the aperture. This assumption is only valid when the pixels are uncorrelated 
and can be treated as independent experiments. Sub-pixel shifts used to register the images, 
astrometric re-projection and PSF matching introduce correlations between pixels. Previous 
work has shown that assuming Poisson statistics in the pr esence of correlated pixels signifi- 



cantly underestimates the u ncertainties in the background (ILabbe et al.ll2003l ; iGawiser et al. 



2006al : iQuadri et al.ll2007bl ). Proper characterization of the noise properties of our images is 
of great importance since the quality of % 2 fits during photometric redshift estimations relies 
on the photometric uncertainties. 



We use the "empty aperture" method (ILabbe et al.l 120031 ) to empirically determine the 
noise properties of our NIR sub-fields and optical images. Briefly, we randomly place a 
large number (5000) of apertures of area n P i X on our images, reject all apertures falling on 
sources, and measure the flux in the remaining ones. A histogram of the measured fluxes is 
constructed and we perform a Gaussian fit to it. The rms of the background in an aperture of 
area n pix is given by the a of the Gaussian fit. We repeat the process six times for each image 
and take the average a. The method is illustrated by Figure [21 which shows the histogram 
of enclosed fluxes for apertures of two different sizes. It is clear that the background noise is 
well described by a Gaussian distribution. Tables [2] and [3] give the 5a point-source limiting 
depth of our images for an aperture size of 1.4xFWHM (see Section 5.2). 

We used the same method to measure the dependence of noise with aperture size on the 
final K-band mosaics (built as explained in Section 4.5). We measured the rms on a series 
of increasing apertures with diameters between 0.5" and 4.0" (see Figure [3]) and fitted the 
dependence using the following parametrization, 



a n = aain; 



(10) 



where u\ is the pixel-to-pixel rms, n pix is the number of pixels in the aperture and a, 
(3 are free parameters. In the Poissonian case we expect a = 1 and (3 = 0.5. We obtain 
a = 1.84 in both fields and (3 = 0.67 in 1030 and 0.71 in 1255 respectively. Our result 
is consistent with previous works showing that assuming Poisson noise in the presence of 
correlated pixels dramatically underestimates photometric uncertainties. 

We built an rms map for each K-band sub-field by taking the inverse square root of the 
normalized exposure map and scaling it to the rms measured in the 1.4xFWHM diameter 
aperture (the choice of this aperture is discussed in Section 5.2). We properly combined 
these maps in order to build a mosaic rms map for the whole field and use it to obtain the 
photometric uncertainty in a given aperture when performing photometry on our sources. 



4.4. Image Mosaicing 

After the astrometric re-projection of the K-band sub-fields we obtained images with 
the same size and physical coordinate system as the MUSYC optical images, so we just 
needed to combine the images of each sub-field in a convenient way in order to create a final 
K-band mosaic of each field. No registration of the frames was necessary. We calculated a 
weight for each pointing given by 



3 (flux scale j ■ rmsj) 2 

The difference of this weight with the one used in Section 4.1, is that it does not take into 
account the FWHM since the sub-fields have been PSF matched, and the rms is no longer 
the pixel-to-pixel rms, but it corresponds to the rms of the background in a 1.4 x FWHM 
diameter aperture estimated as described above. The flux scaling factors correspond to the 
ones used to take all the sub-fields to a common zero-point (see Section 4.3). The normalized 
exposure map of each sub-field is scaled to the corresponding weight in order to create a final 
weighting map. Finally, a weighted average of the nine sub-field images is taken at every 
pixel position in order to create the K-band mosaic. Figure H] shows the final K-band mosaic 
of 1030 and 1255. 
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5. Source Detection, Photometry and K-selected Catalogs 

5.1. Detection 

Source detecti on was performed on the rms normalized K-band mosaic of each field using 



SExtractor v2.5.0 (IBertin fc Arnoutslll996l ). The detection image was built by dividing each 
K-band mosaic by its normalized rms map and then masking the noisy edges of the image 
(product of the dithering) as well as the regions around bright saturated stars. After masking, 
the K-band images cover an area of 830 arcmin 2 in 1030 and 828 arcmin 2 in 1255. 

In order to find an optimal set of detection parameters (DETECT_THRSH, the detection 
threshold in pixel-to-pixel rms units, and DETECT_MINAREA, the number of contiguous 
pixels that must meet the threshold), we estimated the number of spurious sources in our 
catalogs by running SExtractor with the same set of parameters on the negative of the 
detection image, which is a robust estimation given the gaussian noise properties of the 
background. We looked for a set of parameters that maximizes the number of detected 
sources, while keeping the fraction of spurious sources at K < 20 under 1%. We found 
that filtering the images with a PSF sized kernel lowered the number of detected sources 
(keeping our spurious fraction constrained) so we did not apply any kind of filtering during 
the detection. The optimal detection parameters found were DETECT_MINAREA=1 and 
DETECT_THRSH=3.8 (3.6) for 1030 (1255). 



5.2. Photometry and Color Estimation 

Figure [5] shows the normalized Signal-to- No ise (S/N) ratio as a fun ction of aperture size 



for a point-like source in the 1030 field. Following lGawiser et al.l (j2006al ) the S/N ratio is built 
using the growth curve of the Moffat profile to which the PSF was matched and the noise 
dependence with aperture size (u n ( n P ix))- For 1030 and 1255, the S/N ratio is maximized 
by an aperture of ~ l.lxFWHM. The choice of the aperture used to calculate the colors 
of our sources is a delicate issue. Choosing a small aperture of l.lxFWHM maximizes the 
S/N, but makes the photometry highly dependent on the errors in the PSF matching and 
astrometric re-projection of the K-band images. We decided to take a more conservative 
approach and use a "color aperture" of 1.4xFWHM (given by the half-light radius of the 
PSF) for which the S/N ratio is still above 95% of its maximum. 

We performed photometry using SExtractor in dual image mode. In this mode, the 
detection is performed in the K-band detection image and photometric measurements are 
obtained on the images in all bands. We measured the total K-band flux of our sources 
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using a Kron auto-scaling elliptical aperture (SExtractor's FLUX_AUTO). It has been shown 
that Kron apertures mi ss a small fraction (<10%) of the total flux of extended objects 



( iBertin fc Arnoutsl Il996l ). We accounted for this loss by applying an aperture correction 
computed using the growth curve of the PSF and correcting by the fraction of the flux 
falling outside the characteristic radius of the Kron ellipse given by r kron = \fab (with a, b 
the major and minor semi-axis). Corrections are smaller than 0.1 mag for K<20 sources. 

Finally colors were estimated using the flux measured in a fixed 1.4xFWHM diameter 
circular "color aperture" in all bands. Colors, together with the total K-band flux allow the 
determination of total fluxes in all bands. 



5.3. K-selected Catalogs 

K-band selected catalogs are publicly available at \http://www. astro. yale. edu/musyc^ 
All fluxes and uncertainties are given in units of 0.363 /iJy so the zero-point for AB mag- 
nitudes is 25 mag. Photometric uncertainties given in the catalogs correspond to the back- 
ground rms in the corresponding aperture (o~bkg)- Total photometric uncertainties should be 
calculated as 



a phot = yJ(c bkg ) 2 + FLUX (12) 

Version 1.0 of the catalogs has the following format: 
Column 1: SExtractor ID number 
Columns 2-3: x and y centroid (pixels) 
Columns 4-5: a and 5 (J2000.0, decimal degrees) 
Column 6: Internal field code: 1=1030 2=1255 
Columns 7-20: Flux density and error in color aperture (UBVRIzK) 
Column 21-22: K-band aperture-corrected total flux density and error 
Column 23-29: Exposure time weight, normalized to maximum for the field (UBVRIzK) 
Column 30: Color aperture diameter (arcsec) 
Column 31: AUTO aperture diameter 2\fab (arcsec) 
Column 32: Semi-major axis of AUTO aperture (pixels) 
Column 33: Semi-minor axis of AUTO aperture (pixels) 
Column 34: AUTO aperture Position Angle (decimal degrees) 
Column 35: AUTO aperture correction factor 
Column 36: SExtractor star/galaxy class 
Column 37: Object was originally blended with another 
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Column 38: Object's neighbors may significantly bias AUTO photometry 
Column 39: SExtractor flag for the K-band 



6. Star-Galaxy Classification 



We separated sta rs from galaxies in our catalogs using the two color criteria introduced 
by iDaddi et al.l (120041 ) in which stars populate a region in the (B — z)ab versus (z — K) AB 
plane (hereafter: BzK plane) given by [z — K)ab < 0.3(5 — z)ar ~ 0-5. This me thod was 
calibrated using the spectroscopically complete K20 survey ( ICimatti et al.l l2002bl ) and has 
been shown to be extremely ro bust, showing a ve ry high efficiency (>90%) at selecting stars 
and very little contamination (IDaddi et al.l 120041 ) . Since the BzK method was developed by 
Daddi using a slighltly different set of filters than the ones used in this work (the biggest 
difference is the use of a Bessel B-band as opposed to the Johnson B-band used here) we used 



the (B -z) and (z-K) colors of 35 stars in the K20 GOODS-South field given by IDaddi et al. 



( 120041 ) and compared them to the colors of the same star s obtained from our B, z and K-band 
photometry in the ECDF-S field (ITaylor et aLllin prep.1 ). We apply the computed offsets of 



-0.04 mag in [z — K) and 0.56 mag in (B — z) to our colors before plotting them in the BzK 
diagram. 

This method works better on our data than the star-galaxy classification based on the 
SExtractor CLASS_STAR parameter. The lower panel of Figure [6] shows the BzK diagram 
for sources brighter than K = 20 in the 1 255 field . We have overploted stars of different 
spectral types (from 05 to M6) taken from IPicklesI ( 1998 ). all of which fall in the selection 
region. From the lower panel in Figure [6] we can conclude that the CLASS_STAR parameter 
works correctly on our data up to K ~ 16. At fainter magnitudes the method loses its ability 
to separate stars from galaxies. 

After star-galaxy separation we detected 1344 stars and 8015 galaxies in 1030 and 1547 
stars and 8115 galaxies in 1255 to a magnitude limit of K = 20. The number of stars 
detected is in roug h agreement with models of stellar population synthesis in the Galaxy 
(IRobin et al.l 120031 ) which predict a number of ~ 1100 stars brighter than K = 20 for the 
areas and galactic coordinates of our fields. 



7. Number Counts 



We computed differential K-band number counts in 0.5 magnitude bins for all galaxies 
in both fields separately. The results are presented in Figure [7] together with data from 
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the UKIRT Infrar ed Deep Sk y Surv ey (UKIDSS) Ultra Deep Survey Early Dat a Release 
(UDS EDR) from lLane et all (l2007h. the Deep3a-F field from iKong et all hood ) and the 
MUSYC Deep NIR Survey (IQuadri et al.l l2007bl ) . Error bars assume Poisson statistics and 
underestimate the real errors because cosmic variance is not taken into account. The 1030 
and 1255 fields show an excellent agreement (within la) in their number counts in the 
16 < K < 20 ra n ge. W e do not observe the excess in the galaxy number counts observed 
by IQuadri et al.1 (j2007bl ) in the 1030 field, showing how larger areas (the MUSYC Wide 
fields are nine times larger than the MUSYC Deep fields) help to overcome the problem of 
cosmic variance. We do, however, observe an excess of passive galaxies at z ~ 2 in this 
field as reported in Section 9.1. The flux distribution of galaxies in our sample shows an 
excellent agreement with previous works. From the comparison with the UKIDDS UDS 
and the MUSYC Deep NIR Survey, both of which are ~ 1 mag deeper than our survey, we 
conclude that our sample of K-selected galaxies ~80% complete in the faintest magnitude 
bin (19.5 < K < 20). 



8. The BzK selection technique 



Daddi et al.l (120041 ) introduced a two-color selection criteria based on B, z and K-band 
photometry to select star-forming and passively evolving galaxies in the 1.4 < z < 2.5 range. 
This method is known as the BzK selection technique, and presents a series of advantages 
like being almost unbiased against highly reddened galaxies and being able to separate 
star-forming from passive galaxies. Daddi used the 94% spectroscopically complete sample 
of 311 K < 20 galaxies from th e K20 survey which covered 32 arcmin 2 in the GOODS- 



South field (jCimatti et al.ll2002bl ) to calibrate the method. Star-forming galaxies at z > 1.4 
(sBzK) occupy a well defined region in the BzK plane to the left of the solid line shown in 
Figure [8] defined by BzK = (z — K)ab — (B — z)ab > —0.2. Old and passively evolving 
galaxies at z > 1.4 (pBzK) are located on the upper right corner of the BzK plane, to 
the right of the solid line and above the dashed line in Figure [S] in a region defined by 
BzK < —0.2 fl (z — K) AB > 2.5. As stated in Section 6, stars are clearly separated from 
galaxies (especially from those with z > 1.4), and are confined to the region under the 
dashed-dotted line in Figure [BJ defined by (z — K)ab < 0.3(5 — z)ab ~ 0.5. Daddi showed 
that the method is highly efficient at selecting z~2 galaxies and presents little contamination 
from low-z interlopers. The reddening vector in the BzK plane is approximately parallel to 
the sBzK selection criteria, which ensures that the method is not biased against heavily 
reddened dusty galaxies. 
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9. Results 



9.1. BzK Galaxies in the MUSYC Wide NIR Survey 

We have used our data in the 1030 and 1255 fields together with the MUSYC Wide 
data in the ECDF-S to select BzK galaxies over a total area of 0.71 deg 2 to a limiting depth 
of K < 20. A K-band selected catalog was produced for th e ECDF-S us i ng very s imilar 



techniques to the ones used here, and will be presented by iTaylor et al.1 (jin prep J ). By 



combining the three MU SYC Wide N I R Su rvey fields, this is the largest area ever surveyed 
for BzK galaxies. Since iDaddi et al.1 (120041 ) used a slightly different set of filters than the 
ones used here we found the need to apply a correction to our color photometry in order 
to be able to use the BzK selection criteria consistently. We have used the (B — z) and 
(z — K) colors of stars i n the K20 /GOODS field (which is completely covered by our ECDF- 
S imaging) provided by IDaddi et al.1 (120041 ) to compute offsets that we then applied to our 
measured colors. After correcting our photometry we applied the BzK criteria to our data 
to produce a large sample of z ~ 2 star-forming and passive galaxies. 

Figure M presents the positions of all K < 20 sources in the MUSYC Wide NIR Survey 
in the BzK plane. 291 objects show no detections to the la level in the B and z bands, 
appearing as lower limits in (z — K) and having an undetermined (B — z) color in Figure 
|8j They correspond to 1.0% of the detected sources. This fraction agrees with our expected 
1% of spurious sources (see Section 5.1) so we assume most of these objects correspond 
to spurious detections and eliminate them from our analysis. Galaxies falling in the sBzK 
region and showing lower limits in (B — z) as well as those falling in the pBzK region with 
lower limits in (z — K) cannot be unambiguously classified. There are 73 of these unclassified 
objects and for simplicity we leave them outside any further analysis. Including these objects 
does not affect any of our following results in a significant manner. Out of the 24399 detected 
galaxies, we have unambiguously selected a sample of 3261 BzK galaxies (2502 sBzK and 
759 pBzK) at z ~ 2. This is the largest existing sample of BzK galaxies to this depth. 

Over the whole surveyed area sBzK galaxies have a sky density of 0.98±0.06 arcmin -2 , 
where we have used the relative deviation between our fields to include cosmic variance 
when estima ting the erro r . Thi s value is lower than the previous value of 1.2±0.05 arcmin -2 



obtained by iKong et a 



(2006) and consistent with the 0.94±0.17 arcmin 2 value for the 



K20/GOODS sample (IDaddi et al.l 120041 ) (both values showing only Poisson errors). For 
pBzK galaxies we obtain a sky density of 0.3±0.1 arcmin -2 (again including cosmic variance 
in the error), consistent with the values obtained by Kong (0.38±0.03 arcmin -2 ) and Daddi 
(0.22 ± 0.08 armin -2 ). Due to the considerably larger area surveyed in this work compared 
to these previous studies, and the fact that both sBzK and pBzK galaxies are strongly 
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clustered as we show below, our values for the sky density of these objects are much less 
affected by cosmic variance and are therefore more representative of the true sky density of 
these objects. 

Table [4] shows the number of sBzK and pBzK galaxies in our sample and their sky 
densities for each of our three fields. T he GOODS-S field seems to be underdense in 
DRG at z > 2 5 JMarchesini et all l2007h and in optically br ight AGN at high redshifts 
(iDwellv fc Pagell2006h . furthermore, Ivan Dokkum et al.l (120061 ) found that the sky density 
of massive (M > 1O 11 M ) K-selected galaxies at 2 < z < 3 in the GOODS-S field is a factor 
of 3 lower than that of the 10' x 10' MUSYC Deep 1030 field (the densest field in that work). 
Here we find a similar trend for passive galaxies in the 1.4 < z < 2.5 range, where the sky 
density of pBzK galaxies in the 30' x 30' ECDF-S is a factor of 3 lower than in the 30' x 30' 
1030. The density of passive z ~ 2 galaxies in the ECDF-S field is 60% of the mean value for 
the whole survey, while in the 1030 field it is a factor 1.7 higher. Interestingly, star-forming 
galaxies at z ~ 2 show no sign of a significant underdensity. 

The large number of objects in our sa mple allows us to c learly identify the new branch 
of galaxies in the BzK diagram reported by lLane et al.l (120071 ). This feature runs parallel to 
the stellar sequence in the 2.5 < (B — z) < 5.0 range at [z — K) ~ 1.0. Lane found it to 
be consistent with the track of passively evolving early- type galaxies at z < 1.4. This is the 
first confirmation of the existence of this new branch since it was reported. 



9.2. Redshift Distribution 

Photometric redshifts rely on the presence of continuum spectral features in the galaxy 
SED strong enough to show up in broad-band photometry. The most prominent features 
of this kind are the Lyman break at 912 A and the 4000 A break. At z ~ 1.5 the 4000 A 
break is just red-wards of the z-band and the Lyman break is ~1000A blue-wards from the 
U-band. As we move to higher redshifts we have no strong spectral features being sampled 
by our broad-band photometry (UBVRIzK) until we reach z ~ 2.6 where the Lyman break 
enters the U-band (the 4000 A break does not get sampled by the K-band until z ~ 4). This 
means that using our data alone we are unable to estimate reliable photometric redshifts for 
our sample of BzK galaxies since they inhabit a nearly identical redshift range to that in 
which we cannot sample any spectral continuum features. We would not have this problem 
if we had J and H-band imaging of our fields. 

To overcome this problem we have used the data from the MUSYC Deep NIR Survey, 
which covers a smaller area to a greater depth and includes JHK near-infrared imaging, as 
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well as the UBVRIz optical coverage. All the data in the Deep Survey was taken using the 
same instrument as the data presented here, and was reduced using analogous techniques. 
We have limited the Deep Survey catalog to K < 20 and used the same color corrections 
and selection criteria explained above to select a sample of 514 BzK galaxies (365 sBzK and 
149 pBzK). This sample should be representative of our BzK galaxy sample since there are 
no significant differences in the selection biases involved. 

Photom etric redshifts of the Deep Survey galaxies were calcula ted with the methods 
described by iRudnick et all (120011 ) using linear c ombinations of four Coleman et all (119801 ) 
empirical templates, two starburst templates from lKinney et al. (|l996h and two solar metal- 



licity, dust-free lOMyr and lGyr old single stellar populatio n templates (Bruzua 



20031 ). For more details on the photo-z estimations refer to lQuadri et al. 



(I2007bh . 



fc Chariot 



Figure M shows the redshift distribution N(z) of K < 20 pBzK and sBzK galaxies in the 
Deep Survey. Both distributions are fairly Gaussian and cover the expected redshift range 
for BzK galaxies, roughly 1.4 < z < 2.5. The distribution of sBzK galaxies is clearly broader 
than the distribution for pBzK. The peak at the 1.25 < z < 1.5 bin in the pBzK distribution 
comes from the 1030 field and is not present in any of the other three Deep Survey fields. 
As can be seen from Table H] the 1030 wide field shows an excess of pBzK galaxies of a factor 
of 1.7 with respect to the whole survey. This overdensity, together with the presence of the 
peak in the redshift distribution might indicate the presence of a large scale structure at 
z ~ 1.4 in the 1030 field. 



Similar to lHayashi et all (120071 ) we obtain best Gaussian fits of our redshift distributions, 
but in our case we correct them for the effects introduced by the large errors involved in 
photometric redshift estimations as explained below and use these to recover the correlation 
length from the angular correlation function through the inverse Limber projection (see 
Section 9.5). In that work Hayashi used the spectroscopic redshift distribution of bright 
(K < 20.1) sBzK galaxies to compute correlation lengths (ro) for a much deeper sample 
of sBzK (K < 21.3) under the assumption that N(z) is not d e penden t on K magnitude. 
This latter assumption was proved to be wrong by Quadri et ali (2007b) us ing the MUSYC 
Deep NIR Survey data. By looking at Figure 12 of iQuadri et al. (j2007b ) it is clear that 
the redshift distribution of sBzK is severely broadened when going to fainter magnitudes 
(K < 21) as the BzK technique loses its power to select galaxies in a well determined 
and narrow reds hift range. This bro adening of N(z) for K > 20 BzK galaxies had also 
been reported by iReddy et all (120051 ) using spectroscopic redshifts of sBzK galaxies in the 
GOODS-N field. Assuming a much narrower redshift distribution during the clustering 
analysis should underestimate the real value of ro, and this might affect the level of clustering 
measured by Hayashi et al. for their faint sample of sBzK galaxies. 
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The observed N(z) is strongly affected by errors in the photometric redshifts. In 
the redshift range of interest the Deep Survey photometric redshifts show a dispersion of 
Az/(1 + z) = 0.12 when compared to spectroscopic redshifts drawn from the literature 



(IQuadri et al.ll2007al ). To a first approximation, under the basic assumption that N(z) is a 
normal distribution and that photometric redshift errors are of a Gaussian nature, the net 
effect of the errors is to broaden the intrinsic redshift distribution. The observed distribution 
in this case is the result of the convolution of the intrinsic distribution and the photometric 
error distribution. This is a delicate issue because of the assumptions involved in this ar- 
gument. First, we are assuming that the scatter in the photometric redshift is not redshift 
dependent, which given the narrow range over which BzK galaxies are selected seems like a 
reasonable approximation. Also, uncertainties in photometric redshifts are subject to a series 
of systematics related to the methods used to estimate them (templates used, photometric 
errors in different bands, aliasing of SED features, etc) which can lead to a non-Gaussian 
distribution of errors. Nonetheless, unless a very specific and unlikely combination of error 
distribution and redshift dependence of this distribution occurs, the net effect of photometric 
uncertainties is to broaden the intrinsic redshift distribution considerably and this effect must 
be taken into account in order to avoid significant overestimation of the level of clustering 
of the galaxies in the sample. 

The best Gaussian fits to the observed N(z) are centered at zq = 1.58 ± 0.04 and 
z = 1.78 ± 0.03 and have a width of a Zj0hs = 0.35 ± 0.04 and o Zfihs = 0.45 ± 0.03 for pBzK 
and sBzK galaxies respectively, and are shown as the dotted lines in Figure M Errors in 
the parameters come from Monte-Carlo simulations of 1000 renditions of N(z) where the 
redshift of the galaxies were varied within their uncertainties. We have deconvolved these 
fro m our assumed Gaus sian distribution of photometric errors of width Az, as determined 



by IQuadri et al.l (j2007al ). measured at the center of the distribution for each population, to 
obtain our estimates of the intrinsic redshift distribution N corr (z). The intrinsic distribu- 
tions are shown as dashed lines in Figure [HI normalized to the total number of objects in 
each sample, and have widths of a z = 0.17 ± 0.06 and a z = 0.31 ± 0.04 for pBzK and sBzK 
galaxies respe c tively . The latter value is comfortingly consistent with the 0.35 obtained by 



Hayashi et al.l (120071 ) for a sample of 81 sBzK galaxies with spectroscopic redshifts, support- 



ing the validity of the above assumptions. 



9.3. Stellar Masses, Reddening and Star Formation Rates 



Besides calibrating the BzK method, Daddi used iBruzual &: Charlotl (120031 ) models to 
perform SED fitting of the K20/GOODS galaxies with known spectroscopic redshifts to 
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calibrate relations that estimate the stellar mass of BzK galaxies and the reddening of sBzK 
galaxies based on the B,z and K band photometry alone. He also gave a recipe to estimate 
the ongoing S FR o f sBzK using the B flux (rest-frame UV flux) based on calibrations by 



Madau et al.l (119981 ) . F or details on t he ca libration of these methods and the uncertainties 



involved in them refer to lDaddi et al.l (120041 ) . We have used these relations to calculate stellar 
masses, reddening (E(B-V)), and SFR of the galaxies in our sample. It is worth noting that, 
since we do not have redshift information for individual galax ies, these results sh ould be 
interpreted in a statistical sense, especially the SFR estimates. iDaddi et al.l (120041 ) showed 
that mass estimates have uncertainties of ~60% for single objects and E(B-V) estimates 
show a residual rms of about 0.06 when compared to the values obtained from the multi- 
band SED fitting of objects with known spectroscopic redshifts. In any case, given the large 
number of galaxies in our sample and the narrow Gaussian shape of the redshift distributions 
presented in the previous section, the median stellar masses, E(B-V) and SFR should be 
highly representative of the properties of typical BzK galaxies at z ~ 2. 

The top panel in Figure [10] shows the distribution of stellar masses of sBzK and pBzK 
galaxies in our sample. It can be seen that sBzK galaxies span a broader range in mass than 
pBzK, extending down to masses of a few 1O 1O M , while the lowest mass pBzK galaxies in 



our sample have ^ 
a median mass of 



10 M . The median mass of sBzK is ~ 1.0 x 1O U M , while pBzK have 
- 1.6 x 10 11 M Q . These results are in excellent agreement with the results 



of iKong et al.l (120061 ) in the Deep3a-F. As noted in that work, the lower limit of the mass 
distribution of pBzK galaxies corresponds to the mass of a K = 20 galaxy lying on top of 
the (z — K) = 2.5 selection limit, and hence has its origin in the construction of the pBzK 
sample. 

The E(B-V) distribution for sBzK galaxies is shown in the central panel of Figure [TTJJ 
Daddi's recipe to estimate reddening uses the (B—z) colo r as a m easure of the UV slope of the 
SED and assumes a gray and self-similar ICalzetti et al. (hoooh extinction law. The median 
reddening of our sBzK galaxies is E(B-V)=0.47. lAdelberger fc Steidell (120001 ) present the 
distribution of extinction at 1600 A for Lyman-break galaxies at z ~ 3. Under the Calzetti 
extinction law assumption it can be seen that pr a ctical ly no galaxies with E(B-V)>0.45 
are selected by the LBG technique. iReddy et al.l (120051 ) shows that for spectroscopically 
confirmed BM/BX galaxies at z > 1 typical E(B-V) values are lower than 0.3. This means 
that at least 55% of the sBzK galaxies in our sample would not be detected by rest-frame 
UV continuum based selection techniques. 

SFR of sBzK were estimated using the B-band flux as the 1500A rest-frame flux (the 
B-band effective wavelength corre sponds to 1845A at z=1.4, 1593 A at z=1.78 and 1265A 
at z=2.5) and the calibration by IMadau et al.l (I1998I ). UV fluxes were dereddened using 
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the estimated E(B-V) values for each galaxy. As stated above, the lack of knowledge of 
the distance modulus of our galaxies means that these SFR estimates are not reliable for 
individual objects, but the median SFR of the sample should correspond to the typical SFR 
of sBzK galaxies. The distribution of SFR is presented in the bottom panel of Figure [TUl 
and it shows a median val ue of ^230 M^yr - 1 . T his value is s imilar but slightly higher 
than previous estimates by iDaddi et al.l (12004 ) and IfCong et al.l ( 120061 ). Given our survey 
comoving volume of 8.7 x 10 6 Mpc 3 (roughly 1.4 < z < 2.5) we compute a number density of 
sBzK galaxies of ~3x 10~ 4 Mpc" 3 , which translates to a SFRD of ^O.OTM^yr- ^Mpc- 3 . The 
latest estimates of the SFRD at z ~ 2 based on rest-frame UV selected samples (IReddy et al. 



20071 ) show a SFRD~0.2 Moyr^Mpc" 3 at z = 2. NIR bright sBzK can account for about 



a 30% of this total SFRD. Considering that UV selection methods miss aproximately half 
sBzK galaxies, including bright (K < 20) star-formig BzK galaxies in the SFRD would 
increase its value by aproximately a 15%. 



9.4. Number Counts of BzK Galaxies 



In Figure [TT] we present differential K-band n umber counts for s BzK and pBzK galaxie s 



in our sample, together with previous estimates by lKong et al.l ( 120061 ) and lLane et al.l (120071 ). 
The number counts of sBzK galaxies are very steep and the number of objects increases 
sharply towards higher magnitudes. On the other hand, passive BzK galaxies exhibit a 
flattening in the number counts towards fainter magnitudes with a knee around K ~ 18.5. 
This behavior has been observed in previous works and is attributed to the fact that because 
of the narrow redshift distribution of pBzK, the number counts actually probe the luminosity 
function of these objects. Star-forming BzK galaxies do not show this behavior. This has 



(Kone et al. 


2006; 


Lane et al. 


2007) 



pBzK in our sample, this difference in the number counts might as well be related t o intr insic 



differences between the luminosity function of both populations. iMarchesini et al.l (120071 ) has 
shown that the faint-end slope of the luminosity function of galaxies with blue rest-frame UV 
color in the 2 < z < 2.5 range is much steeper than for galaxies sho wing a redder rest - frame 



UV color. A qualitative comparison of Figure [TT] and Figure 5 in IMarchesini et al.l (120071 ) 
supports this statement. Further study is necessary to decouple the effects of these differences 
in the redshift distribution and the luminosity function on the shape of the number counts. 

Number coun ts of pBzK galaxies show an excellent agreement with previous results by 
Lane et al.l (l2007f) on the UKIDSS UDS field, and a somewhat worse agreement with the 
Kong et al.l (120061 ) sample which can be attributed to the small number of objects in the 
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Deep3a-F, especially at bright magnitudes. On the other hand number counts of sBzK in 
our sample are in perfect agreement with the sBzK sample on Deep3a-F (a factor of ~3 
more numerous than the pBzK sample in the same field). Both Kong's number counts and 
the ones presented here are systematically lower than the number counts of sBzK galaxies 
in the UKIDSS UDS sample over the whole magnitude range. The origin of this discrepancy 
lies in the fact that in this work, as well as in Kong's, BzK colors were corre cted in order 



to tak e into account the difference in filters used relati ve to the ones used by iDaddi et al. 
( 12004J ) to calibrate the BzK method. lLane et al.l (120071 ) did not apply this correction, and 
hence their BzK selection criteria was not used consistently with the original calibration of 
the method. Not applying this correction causes an effective offset of ~0.5 mag towards 
lower values in the (B — z) colors of their objects (or equivalently an offset of the selection 
criteria limit towards higher values of (B — z)). As can clearly be seen by inspecting Figure [H] 
this offset will produce a significant excess of objects selected as sBzK but will not affect the 
pBzK sample in such a radical manner (explaining our agreement with Lane's pBzK number 
counts). This excess of low-z contaminants in the UKIDSS UDS sBzK sample can account 
for the discrepancy observed in the number counts. 



9.5. Clustering of BzK Galaxies 

9.5.1. The Angular 2-Point Correlation Function 

The angular 2-point correlation function u{6) is a useful statistic that helps us charac- 
terize the clustering properties of galaxies on the sky plane, in terms of the joint probability 
5P = N 2 [l + uj(6)]5Qi5Q2 of finding two galaxies in the infinitesimal solid angles SQ± and 
8Q2 separated by an angular distance 9 given a surface density of obje cts N, with respect to 



what is expected for a Poissonian random distribution (jPeeblea 119801 ). 



Usually, u)(6) is estimated by comparing the count of galaxy pairs with angular sepa- 
ration 9 found in the observed field to the number of similar pairs obtained from a catalog 
of random and independently distributed objects. Several estimators have been designed to 



quant ify this statistic, and in this work we use the estimator proposed by lLandy fc Szalay 



(119931 ). which has the advantage of being unbiased and showing nearly Poisson variance: 



DD{6) _ 2DR{6) _|_ RR{9) 

r,(Q\ n dd ™dr «rr /1Q\ 

w w mm ( u ) 

Here, DD(6) is the observed number of galaxy pairs with separations in the interval 
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[9,9 + 59], RR(9) is the analogous number obtained from a random catalog that completely 
imitates the field angular geometry, and DR(9) is the amount of observed-random cross 
pairs with separations in the cited range. Each of these quantities is normalized by the 
total number of pairs in the sample, which correspond respectively to tim = N d (N d — l)/2, 
n rr = N r (N r — l)/2 and nd r = N r Nd, with Nd and N r the total number of objects in the 
observed and random catalogs. 

We calculated u)(9) in angular distance bins of constant logarithmic width Alog{9) = 
0.15. In order to produce random samples that replicate the survey geometry, we constructed 
angular masks over the K-band images used for galaxy detection, tracing the field borders 
and leaving out regions occupied by bright foreground stars or any image artifacts that 
forbid source detection. We then generated catalogs of 1000 random points located inside 
the angular mask, and calculated Cj(6) repeatedly over 25 different catalogs in order to 
amount to a total random sample ~ 100 times larger than the observed one. We counted 
random, observed and crossed pairs in each field separately, and then combined the resulting 
DD(9),RR(9) and DR(9) for each bin to produce the angular correlation function for all 
sBzK and pBzK galaxies found in the three fields of the MUSYC Wide NIR Survey. Since 
our three fields have practically identical areas and depths, the number of galaxies in each 
field is dominated by large scale structure rather than observational limits. Because of this 
we have taken the approach of calculating a single correlation function for the whole survey 
by counting all pairs in the three fields instead of computing a single to{9) for each field and 
averaging the results. We consider our approach to be more statistically robust than the 
latter since low number of objects in an angular bin can translate into very noisy estimations. 
A full derivati o n of this approach of handling multiple fields to calculate u{9) will be given by 



Francke et al.l (lin prep.1 ). The uncertainty in £o(9) was estimated using a jackknife method, 
in which we divided each field into 25 smaller regions, and repeated the estimations of Cj{9) 
eliminating one sub-area at a time. In this way the standard deviation for each angular 
distance bin was obtaine d. For a deta i led st udy on the robustness of the jackknife method 
refer to the Appendix on lZehavi et al.l (120021 ). 



The angular correlation function is typically modeled by a power law of the form uj{9) = 
A^ 1 ^ 1 . The limited size of our sample does not allow an independent and significant 
measure of both the amplitude and the slope of w(9), so in the rest of the analysis we 
assume a value of 7 = 1.8, consistent with slopes measured in faint and bright gal axy surveys 



2006 



(Zeh avi et al.l 120021) . This also allows a direct comparison with previous work (IKong et al. 



Hayashi et al.l 120071 ). Since we are estimating u{9) in a finite region of the sky, we 



are affected by an uncertai nty in the estimation of the background galaxy surface density 
(|Peebleslll980l ; Ilnfantdll994l ) that is corrected by introducing a negative offset so that a more 
accurate parametrization of u (9) is given by 
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u(6)=A u {0^-C) (14) 



This bias is know as the "integral constraint" and, as shown by Roche et al. ( 19991 ). 
be estimated numerically by 



can 



Egg: ri5 ^ 

which renders a value of C = 1.89. Once we have applied this correction, we fit the 
amplitude through minimization of x 2 - For fitting purposes we only consider separations 
in the 15" < 9 < 0.1° range. The lower limit corresponds to half a virial radius (r2oo) of 
a halo of ~ 1O 13 M at z=2, in order to avoid the contribution of the "one-halo" term (i.e. 
clustering signal coming from galaxies that share DM halos). We experimented by setting 
larger lower limits for the bins to be considered for the fit (up to 9 > 40"), and observed 
variations in the final values of the correlation lengths ro of the order of 0.1a, so we are 
confident of not being affected by effects introduced by galaxies inhabiting the same halos. 
The upper limit is taken in order to avoid spurious signal coming from border effects. 

The measured 2-point angular correlation function of passive and star-forming BzK 
galaxies in the MUSYC Wide NIR sample, together with the best fitting power-laws are 
shown in Figure [T2J It is clear that there is a positive correlation signal for both sBzK and 
pBzK galaxies, and that passive BzK galaxies are more strongly clustered than star-forming 
BzK galaxies at z~2. Most points in the correlation function of both sBzK and pBzK are 
consistent with a pure power-law fit within the 2a uncertainties and all of them are consistent 
to a 3a level. Therefore we have not found any significant evidence of structure in u(9) for 
any of the two populations. 

Table |6] reports the best fitting values for the amplitudes of u(0). For sBzK galaxies 
we obtain, = 3.14 ± 1.12 and for pBzK we get A^ = 8.35 ± 1.55 (both in units of 10" 3 ). 
These values are systematically lower but still in statistical agreement wi th the only previou s 



estimate of the correlation function of bright K < 20 BzK galaxies by iKong et al.l (120061 ). 
They obtained values of A u = 4.95±1.69 and A w = 10.40±2.83 for sBzK and pBzK galaxies 
respectively in the Deep3a-F field. The systematically lower amplitudes of the correlation 
functions together with the lower surface den sities (see section 9 .1) measured for both sBzK 



and pBzK galaxies in this work compared to IKong et al.l (120061 ) indicate that the Deep3a-F 



field, 8 times smaller than the area covered in this work, show a level of clustering above the 
cosmic average, demonstrating once again the importance of large area surveys in order to 
overcome cosmic variance. 
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9.5.2. Correlation Lengths, Bias and Dark Matter Halo Masses 

Under the assumption of a power-law form for u(9), the 2-point spatial correlation 
function £(r) also corresponds to a power-law with index —7. Usually, £(r) is parametrized 
in terms of the correlation length r and the slope 7 as 



£(r) 



(16) 



The parameter ro can be c alculated from the angular correlation function using the 
inverse Limb er transformation (iLimberl Il953l : iPeebled Il980l ). which leads to the relation 
expressed by iKovac et al.l (120071 ) 



A. 



r(3ti) r F(z)i?^(z)JW*)M*)<fc 



r(i; 



Jo 00 Nco rr (z)dz 



(17) 



where is the amplitude of u;(#), Da (z) is the angular diameter distance, g(z) is a cosmol- 
ogy dependent expression given by 



9W 



Ho 
c 



l + zf Jl + Q M z + Q A [(1 + 



1] 



and F(z) accounts for the redshift evolution of £(r) described by F(z) 
where e = —1.2 in the case of fixed clustering in comoving coordinates, e 
t iering is fixed in pro 



(18) 

if the clus- 



ter coordinates or e = 0.8 according to the prediction of linear theory 
( jBrainerd et al.lll995l ). In this work we assume the first case of constant clustering in co- 
moving coordinates. Finally, N corr (z) corresponds to the redshift distribution of the studied 
population, which we have modeled as Gaussian distributions with the parameters reported 
in Table [6] (see section 9.2). 

We considered two sources of error in our r calculations: (1) the uncertainty in the 
measurement of (cr ro> A u ), and (2) the uncertainty in the photometric redshifts used to 
construct the redshift distribution of both sBzK and pBzK galaxies (cr r0)(Tx )- We estimated 
these two errors separately and added them in quadrature to obtain the uncertainty in 
the reported r values. In the first case we used a Monte-Carlo approach, varying the A u 
values within their measured uncertainties using a normal distribution of errors, and then 
computed the dispersion of the resultant r values calculated using a fixed N corr (z) with the 
parameters reported in Table [6j In order to estimate the errors coming from the uncertainty 
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in the redshift distribution, we followed the approach of lHayashi et al.l (120071 ) who showed 
that the Limber transformation is insensitive to changes in z , with the width a z of the 
redshift distribution the dominant source of uncertainty. As above, we have estimated the 
uncertainty coming from a z by varying it within its error (assumed Gaussian), but now fixing 
A w . Table [6] independently reports both uncertainties for the two populations, together with 
the measured correlation lengths. 

For sBzK galaxies we measured r = 8.89±2.03 Mpc, value that is much lower, although 
bearly consistent to a la level, than the r = 12.141U Mpc (8. 51^3 ^ioo^P c ) measured by 



Havashi et al.l (120071 ) for K < 20 sBzK galaxies in the Deep3a-F sample by iKong et al. 
r|2Q06h . This is in agreement with the above evidence showing that the Deep3a-F field is 
denser and more clustered than the median field. We also measured r = 10.82 ± 1.72Mpc 
for pBzK galaxies. This is the first measurement of the correlation length of passive BzK 
galaxies ever reported. 

In the context of the ACDM scenario, galaxies form by the cooling and condensation 
of baryonic gas in the cores of dark matter halos. The galaxy spatial correlation function is 
then associated with the auto-correlation function of DM halos and the spatial distribution of 
galaxies of a certain type is biased in t he same way as tha t of their hosting halos with respect 
to the underlying mass distribution (IMo &: White! Il996l ). In t he context of the ellipsoidal 
collapse model extension of the Press-Schechter formalism by ISheth et al.l (120011) they are 
related to fi rst approximat i on by th e linear bias. Following the methods of iQuadri et al. 



( I2007al ) and iFrancke et al.l (lin prep.l ). we estimated the galaxy effective bias adopting the 
following definition: 



a 



8, gal 



,DM 



(19) 



where erf corresponds to the variance (of galaxies or dark matter) in 8/i^QoMpc radius 
spheres. We then estimated the minimum mass M^jj of DM halos hosting galaxies in our 
sample, given by 



f™ min b{M DH )n{M DH )dM 1 



DH 



} eff 



r^min n{M DH )dM A 



(20) 



DH 



where &(Md#) is the bias parameter for halos of mass Mbh taken from 



( 2001 ). and n^Mmj^dMuH is the halo mass function derived by lSheth &: Tormen Jl999h . We 



Sheth et al. 



assumed the simplest case of one galaxy per halo. All the above expressions were evaluated 
at the center of the redshift distribution for each of the two populations studied, that is 
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Z 



1.58 for pBzK and z Q = 1.78 for sBzK, and obtained values of b, 



eff 



2-93185 fOT 



star- forming BzK, which correspond to a minimum halo mass of M™jj = 4 x 1O 12 M , and 



J eff 



3.27±q;47 for passive BzK galaxies, corresponding to = 9 x 10 12 M Q . 



By comparing the bias of different populations of galaxies at different redshifts we can 
obtain insights on the evolution of these systems. Figure [13] shows tracks for the evolution 
of the bias with redshift together with meas ured bias factor value s for different popu l ations 
drawn from the literature. Previous work by lQuadri et al.l (12007at) andlGawiser et al.l (120071 ) 
have presented similar plots using the "no-merging" model of iFry ( 1996 ). As dark matter 
becomes more clustered with time the bias factor of a biased population decreases. It is 
important to note that we expect mergers to have a more important role in higher density 
regions than in less dense regions where galaxies have a lower probablity of merging. Hence, 
including the effects of merging would translate into a steeper evolution of the bias with 
redshift. Therefore, the tracks calculated using the "no-merging" model only provide an 
upper limit for the bias factor of a given point at lower redshifts. In this work we pres e nt bia s 
evolution tra c ks obtained from the "halo merging model" developed by iBond et al.l (Il99ll ); 
Bower! (119911 ) ; lLacey &: Cold (119931 ) in which at any redshift halos have a given probability of 
merging into higher mass halos, so there is a defined mass distribution of descendant halos. 
For a given halo mass (bias) at a given redshift, the tracks follow the bias of the most likely 
descendant halo population, that is, the mode of the conditional mass distribution function 
evaluated at lower redshifts. 

It can be seen that although passive BzK galaxies in our sample show a higher bias 
than star- forming BzK galaxies, their clustering level is still consistent, and both populations 
inhabit halos in a similar mass range. Both sBzK and pBzK galaxies are consistent with 
being the descendant of the b right LBG population at z ~ 3. Bright L BG have typical stellar 
masses of ~ 10 10 - lO n M JPapovich et al.l l200ll : llwata et alil2005h so, given the ~ 2 Gy 
time elapsed between z = 3 and z = 1.6, they could easily be able to assemble enough mass 
to reach the typical masses of BzK galaxies given a reasonable SFR of 10 — 100M Q yr _1 . High 
redshift LAE populations are not consistent with being the progenitors of BzK galaxies. It 
can also be seen that K bright UV selected galaxies at z ~ 2 (the K < 20.5 BX population) 
are more clustered than BzK galaxies, although the BX sample clustering measurement could 
suffer from cosmic variance (lAdelberger et al.ll2005bl ). In any case, the trend is consistent 
with what is observed at fainter magnitude limits, where faint sBzK galaxies show a lower 
clustering level than the BM/BX sample. On the other hand, as can be seen in Figure UM the 
NIR bright BzK galaxies presented in this work are consistent with being the predecessors of 
massive early-type > 1L* galaxies in the local universe usually present in groups and clusters, 
but are totally inconsistent with being associated with the progenitors of the central galaxies 
of rich clusters. The most likely descendants of K bright sBzK and pBzK galaxies have bias 
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factors of l.St° j and 1.6±g;? at z = 0. 



10. Summary and Conclusions 



In this paper we have presented K-band imaging of two of the three 30' x 30' fields 
that conform the MUSYC Wide NIR Survey. We have given details of our data reduction 
procedures, as well as our source detection algorithms. By combining our data with the 
MUSYC UBVRIz optical imaging, we have constructed K-band selected source catalogs 
that reach a 5a point-source limiting depth of K = 20. Catalogs are publicly available as 
part of the MUSYC Public Data Release, and present spatial and photometric information 
for the 16130 galaxies and 2891 stars brighter than K = 20 detected in the 1030 and 1255 
fields. We have also presented K-band differential number counts for galaxies in the two 
fields, which show an excellent agreement with previous K-selected samples. 

By combining our data with the MUSYC K-selected catalog of the ECDF-S field, we 
have selected a sample of 2502 star-forming and 759 passively evolving BzK galaxies at 
z ~ 2 over a very large area, allowing us to study the spatial distribution of these galaxies 
minimizing the effects of cosmic variance. We reported sky surface densities of 0.98 ± 0.06 
arcmin -2 for sBzK and of 0.30 ± 0.10 for pBzK galaxies. The large area surveyed over 
three non-contiguous fields allows the measurement of realistic uncertainties in the reported 
densities. We have found that the ECDF-S field is underdense in passive galaxies at z ~ 2 
when compared to the mean for the complete survey, in agreement with previous studies on 
the density of DRG and K-selected massive galaxies. We have also confi rmed the existenc e 
of the passively evolving early-type galaxy track at z < 1.4 discovered by lLane et al.l (120071 ). 



Using the methods calibrated by lDaddi et al.l (120041 ). we estimated stellar masses of the 
BzK galaxies in our sample and showed that star-forming BzK galaxies have typical masses 
of ~ 1.0 x 10 11 M Q , while passive BzK objects have ~ 1.6 x 1O 11 M . We also estimated the 
reddening and the SFR in our sample of star-forming BzK objects. sBzK galaxies show a 
median reddening of E(B-V)=0.47, implying that at least 55% of them would be missed by 
UV continum based selection techniques. This is in excelle nt agreement with current studies 



on the overlap of ga l axy p opulations at these redshifts (IReddy et al.l 120051 ; iQuadri et al. 



2007bl ; iGrazian et al.l 120071 ) . Star- forming BzK galaxies also showed large median SFR of 
~230 MQyr" 1 , implying that bright K < 20 sBzK can account for up to 30% of the SFRD 
at z ~ 2. 



We also presented redshift distributions of both sBzK and pBzK, which show that both 
populations are selected over a narrow and well determined redshift range, with the pBzK 
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population showing a narrower distribution than that of sBzK galaxies. We corrected these 
distributions to take into account the broadening introduced by errors in the photometric 
redshift estimation, and used them to deproject the angular correlation function of the two 
populations . 

The 2-point angular correlation functions for pBzK and sBzK galaxies were presented, 
together with their best fits. We used the corrected redshift distributions to deproject the 
spatial correlation function from the angular one and estimated correlation lengths r of 
8.89 ± 2.03 Mpc and 10.82 ± 1.72 Mpc for sBzK and pBzK galaxies respectively, which 
translate into bias factor values of 2.93^0 60 an< ^ 3.27!q47 respectively. By comparing the 
effective bias of the spatial distribution of the two populations with that of DM halos, 
we estimated minimum halo masses of 4 x 1O 12 M and 9 x 10 12 M Q , for sBzK and pBzK 
respectively. 

Finally, we compared the bias of different populations of galaxies at different redshifts, 
and concluded that NIR bright K < 20 sBzK and pBzK galaxies are consistent with being 
the descendants of bright LBG at z ~ 3, and the progenitors of local galaxies with < 1L*. 
This corresponds to the bright end of the luminosity function of present-day galaxies. The 
fact that K bright BzK galaxies evolve into very bright and massive systems is not surprising 
given their high stellar masses and large SFR at z ~ 2. In the future, pushing the limits of 
K-band surveys to fainter magnitudes will allow us to detect large samples of the progenitors 
of current < 1L* galaxies through BzK selection, overlapping these samples with the LAE 
population, opening a window to a more complete study of the ancestors of normal galaxies 
in the local universe. 
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Table 1. MUSYC Wide NIR Fields 



Field 


Equatorial Coord. 


E(B-V) a 




a 8 


mag 


SDSS1030+05 


10:30:27.1 05:24:55 


0.02 


Castl255 


12:55:40.0 01:07:00 


0.02 


E-CDFS 


03:32:29.0 -27:48:47 


0.01 


HSchleeel et all 


Jl998r) 





Table 2. 1030 and 1255 subfields K-band image properties 



Field 


Sub-field 


Zero- 


point 


FWHM a 


5a depth b 








mag 




arcsec 


mag 


SDSS1030+05 


C 


21 


.48 


± 


0.02 





.9 


20.22 




NE 


22 


.01 


± 


0.02 


0. 


.9 


19.87 




E 


21 


.99 


± 


0.02 


0. 


.9 


19.87 




SE 


22 


.11 


± 


0.03 


0. 


.9 


20.00 




S 


21 


.91 


± 


0.03 


0. 


.9 


19.69 




sw 


22 


.02 


± 


0.02 


1. 


.0 


19.73 




w 


22 


.05 


± 


0.02 


1. 


.1 


19.92 




NW 


21 


.98 


± 


0.02 


1. 


.0 


19.68 




N 


22 


.03 


± 


0.02 


1. 


.1 


19.72 


Castl255 


c 


21 


.94 


± 


0.02 





.9 


19.86 




NE 


22 


.44 


± 


0.03 


0. 


.9 


19.84 




E 


22 


.25 


± 


0.03 


1. 


.0 


19.72 




SE 


22 


.38 


± 


0.01 





.9 


20.00 




S 


22 


.40 


± 


0.03 


0. 


.9 


19.77 




sw 


22 


.20 


± 


0.03 


1. 


.0 


19.77 




w 


22 


.27 


± 


0.03 





.9 


19.89 




NW 


22 


.40 


± 


0.02 


0. 


.9 


20.04 




N 


22 


.28 


± 


0.02 


1. 


.0 


20.03 



a Before PSF matching. 

b Vega magnitudes. After PSF matching 
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Table 3. MUSYC Optical Images of 1030 and 1255 



Field 


Filter 


Zero-point 


FWHM a 


5cr depth b 






mag 


arcsec 


mag 


SDSS1030+05 


U 


22.06 


1.1 


25.56 




B 


24.97 


1.0 


25.91 




V 


25.48 


0.9 


25.90 




K 


25.75 


0.9 


25.81 




I 


25.36 


0.9 


25.00 




z 


24.36 


1.0 


23.75 


Castl255 


U 


23.32 


1.2 


25.62 




B 


24.84 


1.3 


25.53 




V 


24.90 


1.1 


25.42 




R 


25.66 


1.2 


25.08 




I 


24.94 


1.1 


24.00 




z 


24.30 


1.0 


22.81 



a Before PSF matching. 

b AB magnitudes. After PSF matching 



Table 4. Number of objects and sky densities per field 



Field 



Galaxies sBzK sBzK arcmin 



pBzK pBzK arcmin 



1030 
1255 
ECDF-S 



8015 
8115 
8269 



758 
911 

833 



0.92±0.03 
1.10±0.04 
0.93±0.03 



441 
165 
153 



0.50±0.02 
0.20±0.02 
0.17±0.01 



Total 



24399 



2502 



0.98±0.06 



759 



0.30±0.10 



a Errors for individual fields correspond to Poisson errors. For the total survey we 
have included cosmic variance in the error estimation. 



Table 5. MUSYC BzK galaxies differential number counts 



K bin center sBzK pBzK 



mag log(N deg 2 mag x ) 



16.75 


0.763 




17.25 


0.939 


1.064 


17.75 


1.541 


1.638 


18.25 


2.275 


2.313 


18.75 


2.860 


2.674 


19.25 


3.300 


2.864 


19.75 


3.615 


2.867 
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Table 6. Clustering Properties of MUSYC BzK Galaxies 



Population 


N aai 




4 o 




^0 






Bias 








10~ 3 














M 


pBzK 
sBzK 


759 
2502 


8.35 ± 1.55 
3.14 ± 1.12 


1.58 
1.78 


0.16 ±0.04 
0.31 ±0.06 


10.82 ± 1.72 
8.89 ± 2.03 


1.35 
1.00 


1.07 
1.77 


o 97 +0.46 
'>- zl -0.47 
9 no+0.59 
z.ao_ 60 


1 x 10 13 
4 x 10 12 



a Amplitude of the angular correlation function with slope parameter 7 = 1.8 
b Parameters of Gaussian fit to the redshift distribution N corr (z). 
c Dispcrsion in ro due to photometric redshift errors. 
d Dispcrsion in ro due to uncertainty 
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1.8 2.0 2.2 2.4 

Wavelength (microns) 



Fig. 1. — Transmission curves of the K' (dashed) and the Ks (solid) filters after multiplying 
by the detector's quantum efficiency and the atmospheric transmission. The transmission 
for the Ks 2MASS filter (dotted) is also shown. 
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Fig. 2. — Histogram of enclosed fluxes measured on ~5000 empty apertures on the 1030 NE 
sub-field, for two aperture sizes of 0.6" and 1.1" in diameter. Best fit Gaussians are also 
shown. 
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Fig. 3. — Measured background rms (dots) as a function of the number of pixels in the 
aperture, and best fit (solid line) for the K-band mosaic of the 1030 field. The dotted lines 
show the formal Poisson scaling dependence with a n a^/n pix (bottom) and a limiting case 
with a n an P i X (top). The vertical dashed line corresponds to the aperture used for color 
photometry (see Section 5.2) 




Fig. 4.— Final K-band ~ 30' x 30' mosaics of the 1030 (left) and 1255 (right) fields 
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Fig. 5. — The solid line shows the Signal-to-Noise for a point-like source as a function of 
aperture size for the 1030 field. The dashed curve corresponds to the enclosed flux at a given 
radius for the field PSF (after PSF matching). The dotted line shows the noise as a function 
of aperture size presented in Figure [31 The vertical dashed-dotted line marks the radius of 
the 1.4xFWHM color aperture used for the photometry. 



-43- 



< 
h- 



as™! 



T3 



0.6 

0.4 



£ 0.2 

_j 
o 

0.0 



1 3 



1 4 



1 5 




: ."4 ,. 



16 17 
K (Vega) 



19 



20 



< 
4? 












4 



/AB 



Fig. 6. — Top: K-band magnitude versus SExtractor's CLASS_STAR parameter for sources 
brighter than K = 20 (dots) in the 1255 field. Objects with CLASS_STAR ~1 are supposed 
to be point-like. The dashed lines show the K < 16 and CLASS_STAR>0.95 region where the 
separation works well. Squares correspond to objects classified as stars using the (z-K)ab < 
0.3(1? — z)ab — 0.5 criteria. Bottom: BzK diagram for sources brighter than K = 20 (dots) 
in the same field. Triangles show sources with K < 16 and CL ASS_STAR>0.95. Stars 
correspond to stars of different spectral types from iPicklesI ( 119981 ). The dashed-dotted line 
shows the (z — K)ab = 0.3(5 — z)ab — 0.5 criteria used to separate stars from galaxies. 
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Fig. 7. — Differential K-ban d number counts f o r galaxies in 1030 an d 1255, together wi th 
data from previous works by iKong et al.l (120061 ); iQuadri et al.l (j2007bl ); lLane et al.l (120071 ) 
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Fig. 8. — BzK plane. Dots show the corrected colors of all K < 20 sources in the three 
fields of the MUSYC Wide NIR Survey. The diagonal solid line delimits the sBzK selection 
region {BzK > —0.2). The horizontal dashed line defines the pBzK selection region ((z — 
K)ab > 2.5). The diagonal dashed-dotted line shows the star-galaxy separation criterion 
((z-K) AB <0.3(B-z) AB -0.5). 
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Fig. 9. — Photometric redshift distribution of K < 20 pBzK and sBzK galaxies in the 
MUSYC Deep NIR Survey. Dotted lines correspond to the best Gaussian fit to the observed 
redshift distribution N obs (z). Dashed lines correspond to the distribution after correcting 
for the broadening due to errors in the photometric redshifts. The parameters reported 
correspond to the corrected distribution N corr (z) 



-47- 




350 : 




E(B-V) 



700 ; 




log(SFR [M yr 1 ]) 



Fig. 10. — a) Distribution of stellar-masses of sBZK (solid) and pBzK (dotted) galaxies, b) 
Histogram of E(B-V) values for star-forming BzK galaxies in the sample, c) Histogram of 
SFR calculated using the dereddened B-band flux of sBzK galaxies. Arrows indicate median 
values of the distributions. 
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Fig. 11. — Differential K-band number counts of K < 20 sBzK (squares) and pB zK (circles) 



i n the MU SYC Wide NIR S urvey, together with data from previous works by iKong et al. 
d2006h and lLane etaD J2007h . 




Fig. 12. — Measured angular correlation function of sBzK (triangles) and pBzK (circles) 
galaxies in the MUSYC Wide NIR Survey. Best power-law fits for sBzK (dotted) and 
pBzK (dashed) galaxies are also shown. The vertical dotted-dashed lines corresponds to the 
15" < 9 < 0.1° range used for the fitting. 
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Fig. 13. — Tracks of bias redshift evolution calculated using the "halo merging" model. 
Large filled stars show bias of K < 20 sBzK an d pBzK calculated i n this work. Also 
presented are the bias of: faint sBzK g alaxies fromlHavashi et al.l (120071 ) (small filled star); 
BM and BX galaxies w i th 1Z < 25.2 from lAdelberger et all (j2005al ) an d K < 20 5 BX g alaxies 
from lAdelberger et al.l (l2005bl ) (open triangles); LBG at z ~ 3 from lLee et al.l (120061) (filled 



triangles), K < 21 DRG and K-selected galaxies at 2.0 < z < 3.5 fromlQuadri et al.l (120 073) 
(squar e s), LAE at z=3.1, z=4.5 and z=4.86 (open stars, Gawiser et al. ( 2007 ^ Kovac et al. 



' — ■ — ~ / i -. - 1 \ ~ jr 7 . — — . . ^ M ^^ j -- M .-^j^^ — — 1 

( 120071 ) ; IOuchi et al.l (120031 ). res pectively), local SPS S galaxies from lZehavi et al.l (120051 ) (open 
circles), and rich clusters from iBahcall et al.l (120031 ) (filled circle). 



- 50 - 



m 

D 
CD 




1.0 1.5 
Redshift 



Fig. 14. — Tracks of bias redshift evolution of K < 20 pBzK and sBzK galaxies. Dotted 
lines show the bias of the most likely descendant halo population. Dashed lines mark the 
bias of halos whose likehood is reduced by a factor exp(— 1/2) with respect to the maximum 
likehood, in the case of a gaussian distribution of descendants this would correspond to ±lcx. 
The probability distribution of descendant halo masses is clearly non-symetric. 



